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Classical and quantum perturbations can be described in terms of marginal distribution functions in the frame- 
work of tomographic cosmology. In particular, the so called Radon transformation and the mode-parametric 
quantum oscillator description can give rise to links between quantum and classical regimes. The approach re- 
sults a natural scheme to discuss the transition from the quantum to the classical perturbations and then it could 
be a workable scheme to connect primordial fluctuations with the today observed large scale structure. 
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I. INTRODUCTION 

Addressing the problem of cosmic evolution, starting from the inflationary mechanism HI Hj] up to the observed scenario, as 
well as the problem of producing classical perturbations starting from the quantum ones, have been the subject of several studies 
during the last decades . In particular, the Friedmann-Lemaitre-Robertson-Walker (FLRW) homogeneous and isotropic 

metric, described by the expansion factor a(t), is usually considered the standard background on which formulate the classical 
and quantum theory of cosmological perturbations. The goal is to provide a suitable description which, starting from quantum 
regime can lead to the observed large-scale structures such as galaxies, galaxy clusters, super-clusters and voids. 

On the other hand, the task of quantum cosmology is to achieve a quantum description of the cosmic initial state which, by the 
following evolution could determine, in principle, all the features of the latest classical epochs which we today observe by the 
astrophysical measurements. The initial quantum state is usually associated with the so-called "Wave Function of the Universe", 
in analogy with standard quantum mechanics Jfl, or with its density matrix J3, HI]. A useful representation to consider the 
quantum-to-classical transition for the quantum state of the universe is the Wigner function W(q,p) (see, e.g., (Hill]). 

Such a Wigner function is the analogue of the classical probability distributions on phase-space f(p, q) which are, in general, 
non-negatively defined. On the other hand, due to the uncertainty relations, the Wigner function can assume negative values in 
some domain of phase-space. In view of this fact, it is not a fair probability distribution and it is called a quasi-distribution. 

Besides, quantum cosmology in its development has always adopted new results from quantum mechanics and quantum field 
theory. Recently, a new tomographic probability representation of quantum states has been found fioll . In this picture, the fair 
nonnegative probability distribution, containing a complete information on the states, is used instead of the wave function or 
the density matrix. Though the new probability representation is essentially equivalent to all the other available representations 
adopted in quantum mechanics fiUl . it has its own merits. An important property of this representation is that the quantum 
state is associated to the same tomographic probability density W(X, /i, v), connected with the Wigner function by the standard 
integral Radon transform [12]. Besides, the classical state is described using the Radon component of the standard classical 
probability density f(p, q) IU3I1 . Thus, in the tomographic probability representation, both classical and quantum states can be 
described by the same non-negative probability densities W (X, fj,,v): this fact makes easier to consider the classical-to-quantum 
mutual relations and transitions. 

Some aspects of quantum cosmology have already been studied in the framework of the tomographic probability approach 
giving a picture coherent with both inflationary and standard Friedmann cosmology ffl [l5l [Trjl ITtL [l8i [l9ll . Since the issue of 
quantum-to-classical transition plays a main important role in the theory of cosmological perturbations, the aim of this paper 
is to study the perturbations using the tomographic probability densities for the description of cosmological "states" (both 
classical and quantum ones). Such "states" are configurations of the universe which can be suitably defined by assigning a 
set of cosmological observables. As shown in Jit], it is possible to find the same equations for perturbations in the cases of 
hydrodynamic matter, of scalar field and of alternative theories of gravity as f(R). The perturbation evolution can be described 
by the Hamiltonian formalism corresponding to small vibrations on a classical background. The field theory can be considered 
in the tomographic picture as shown in [20]. The Hamiltonian for small vibrations corresponds to interacting oscillators with 
time-dependent frequencies and mutual coupling constants. In quantum mechanics, the one-dimensional parametric oscillator 
has been studied by Husimi [21] and new integrals of motion, linear in field quadratures of the oscillator have been found in 
I22I l23ll . A comprehensive approach to the system of interacting parametric oscillators has been presented in J24ll . The theory 
of such oscillators with dissipation can be developed in connection with the non-stationary Casimir effect studied in l29tl . 

The advantage of the tomographic approach for this particular Hamiltonian is related to the fact that, in both quantum and 
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classical domains, the propagators providing the tomograms of states at time t, in terms of the states tomograms at some earlier 
time, are identical. For these particular propagators, the difference between the quantum and classical universe descriptions is 
connected to the choice of the initial state tomogram only. 

The layout of the paper is the following. In Sect. II, the main points of the field theory for cosmological perturbations 
is sketched. The probability representation of cosmological perturbations is discussed in Sect.III. A general scheme for ap- 
proaching classical and quantum perturbations is given in Sect IV. Sect.V is devoted to a detailed description of the one-mode 
parametric quantum oscillator, while the multi-mode case is considered in Sect. VI in view of the classical-to-quantum relations. 
Conclusions are drawn in Sec. VII. 



II. THE FIELD THEORY OF PERTURBATIONS 

In JH, a complete gauge invariant theory for the cosmological perturbations is discussed. The cases of hydrodynamic and 
scalar fields are studied and the formulation in terms of theories alternative to general relativity is considered. 

Gauge invariance is a crucial requirement for quantization, because it allows to properly quantize the physical degrees of 
freedom. 

According to the inflationary paradigm, the origin of the cosmological perturbations, which lead to the hierarchy of structures 
which form the observed universe, can be retraced back to the quantum fluctuations of the very early universe. 

The quantization of cosmological perturbations allows, in principle, to study the initial state of the universe, at least from 
a theoretical point of view. In order to relate the early quantum states of the universe to the cosmological observations, it is 
necessary to develop a model for the evolution of the quantum fluctuations into classical perturbations. 

A unifying picture of these different models for perturbations has been proposed in order to quantize them. 

To this aim, a total action for gravity and matter is 



S = 



1 
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if we are considering General Relativity. On the other hand, actions of the forms 
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s 



-gd x 



F($)R+-gi»'$. lt <fi. v -W($) 



+ / C m (g)V^gd 4 x , 



(3) 



where f(R), F(&), W($>) are generic functions of the Ricci scalar R and a scalar field respectively can be adopted in view 
to face several cosmological problems ranging from inflation to dark energy (see e.g. ll25l I26ll27ll28ll for reviews). 

Considering the simplest case (General Relativity), the action for the perturbations can be derived from (Q~|) by writing first the 
metric (in the conformal time rf) in the ADM formulation 



where the lapse function is 



the shift functions are 



the spatial metric is 



and its inverse is 



ds 2 = (N 2 -NiN^drj 2 - 2Midx'dr] - 7yc&W , 



7 y = a 2 (l - 2ip)6ij + 2a 2 E il3 , 



(4) 
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where the convention of summation has been adopted on the repeated lower indices. 

The fields <j>, ip, E and B, introduced in the previous equations, form the tensor Sg^v which represents the scalar perturbations 
of the FLRW spacetime metric. Each of these perturbations terms are gauge dependent, but as shown in [3], linear combinations 
of these functions and of their derivatives can be taken into account to construct gauge invariant objects. 

Writing the gravitational part of the action in the ADM form and and expressing it in terms of the perturbed metric, we arrive 
at the equation for a gauge invariant combination of the fields <fi, ip, E and B. The perturbations have been studied in three cases, 
i.e. hydrodynamical, scalar field and /(i?)-gravity. In all these three cases, the action for the perturbations, described in terms 
of gauge invariant fields (omitting total derivatives), takes the form 
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Sperturb. = ~ j \ V - C^UiU j + — V ) y/j d X, (9) 

where z is a time-dependent function. Eq. (O is the action of a scalar field v with a time-dependent mass m 2 (77) = — . 
Turning to the Hamiltonian formulation, the conjugate momentum for the field v is 

t(?7,x) = — = 1/(77, x), (10) 



and the resulting Hamiltonian is 



where c s is the sound speed and 



H = \J (» 2 + c h"«,i»j-y» 2 J V7<ft«, (ii) 



Hc s 

The transition to the quantum formulation is obtained once the variables v and ir are replaced with the operators v and n, 
satisfying the following commutation relations 

[v(rj, x), v(n, x')] = [7t(r), x),n(f], x')] = 0, [v(rj, x), -k(t], x')] = itf(x - x') (13) 

where the delta function 5(x — x') is normalized by requiring 

y/^5(x~x')d 3 x = 1. (14) 



III. THE PROBABILITY REPRESENTATION FOR COSMOLOGICAL PERTURBATIONS 

Instead of going into the usual procedure of canonical quantization of the above system, we want to indicate a different 
approach to quantize the perturbation field in terms of probability distribution functions in the way proposed in [ 10]. 

In this formulation of quantum mechanics, marginal distribution functions with classical-like evolution replace the wave 
functions. The main advantage of this formulation is that a quantum theory can be entirely expressed in terms of observable 
functions, which are comparable with their classical counterparts. 

This approach has been already applied in quantum cosmology with the purpose to study the evolution of a quantum universe 
into a classical one and to obtain all the information of the initial quantum cosmological stages from the today observations. 

Recently this approach has been extended to quantum field theory luOll . We shall briefly recall the results here. 

Let us consider the quantum Hamiltonian for a scalar field in a (d + 1) spacetime 

d d x, (15) 
and the combination 



H = 



1 ^ 

— 7T 

2 



~X>0(aO) 2 + tf(0(*)) 
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= [i(x)tp(x) + v(x)fr(x), 



(16) 
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introducing the quantum characteristic function 1 

X(H X )) = ^ ex P J d d xk(x)<f>(x) 
we can define the marginal distribution functional 



(17) 



w($(x),m(x),//(x)) = / me-^ ddxk ^^x{l ' 
which satisfies the following evolution equation ll20ll 
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where Af(x) = f(x + Ax) — f(x) and the operator (— <5/<5$(x)) _1 is defined by 



(19) 
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Vke ->! d d xk(x)<S>[x) _ j vk 1 e -i f d d xk(x)Mx) 

k(x) 



and the dot represents the time derivative. The above formulas allow to develop a theory of cosmological perturbations at 
classical and quantum levels. 



IV. CLASSICAL AND QUANTUM COSMOLOGICAL PERTURBATIONS 

In view of Hamiltonian ( fTTT > and equation ( fT9l ) the evolution of cosmological quantum perturbations can be described in terms 
of a marginal distribution function which satisfies the equation 



W (v(x),^(x),v(x),t) 



<Tx 



s s 

ov(x) OV(X) 



8v(x) ) 8\i{x) 



— (ivjx) 6 
z \ 8fj,(x) 



W(v(x),/j,(x),v(x),t) . 



(20) 



Due to the fact that the potential is quadratic, classical perturbations follow a similar evolution equation. The only difference is 
in the initial conditions which are restricted by the Heisenberg uncertainty principle for the quantum perturbations. A relevant 
expression is also the evolution equation for the Fourier transform of the tomogram 



X(k,li,v,t) = I dXe lkx W(X,fi,is,t) 



(21) 



Note that the meaning of the function Q(x) is different with respect to the scalar field mentioned in the action {3} ■ 
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It is 



d d x 



1 d 

— V\ x ) , x — 2ihk(x)v(x)lS. 
m ov{x) 



ik(x) Sfi(x) 



-2—v{x) 



z 5p,{x) 



X(k(x),/j,(x),v(x),t) . 



(22) 



V. THE ONE-MODE-PARAMETRIC QUANTUM OSCILLATOR 



In the previous section, the Hamiltonian describing the perturbations has been constructed and presented in the form of a 
sum of Hamiltonians of oscillators with time-dependent frequencies. In this section, we are going to consider in detail the one- 
mode evolution for the quantum parametric oscillator both in the Schrodinger representation and in the tomographic probability 
representation. Let us use, for the one-dimensional parametric oscillator, dimensionless units, i.e. with the Planck constant 
h = 1, the "mass" of the oscillator m = 1, and the time-dependent frequency, at the characteristic initial time to, equal to unity, 
i.e. LLi(to) = 1. The Hamiltonian of the parametric oscillator in these units reads 



H 



2 2 

This system has two integrals of motion, linear in position and momentum 1124 13( 



(23) 



A(t) = -j=(e(t)p-i(t)q) 



(24) 



A\t) = --j= (s*(t)p-i*(t)q). 



(25) 



In Eq.(l24b. the complex function of time e(t) (where t can be the conformal time) obeys the classical equation of motion for the 
oscillator 



e(t) +uj 2 (t)e(t) =0. 

The initial conditions 

e(to) = l, i(t )=i, 
provide the commutation relations of the integrals of motion (l24l and (|25T > 

[i(t),it(t)] = l. 



(26) 



(27) 



(28) 



For the initial time to (we shall use the initial time to = 0) the integrals of motion coincide with the standard creation and 
annihilation operators 



A(t ) =a, 



i t (to) = a t 



(29) 



where 



a = -j= (q + ip) 



V2 



(q - ip) 



Here q and p are the position and momentum operators, respectively. The Schrodinger equation 

. . ld 2 lj>(x,t) UJ 2 (t)x 2 . 

has the solution tpo(x, t) which corresponds to the initial vacuum state, obeying the equation 



(30) 



(31) 



i(*)Vo(M) = 0, 



(32) 
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with the initial condition 



ipo(x,t ) = -^=e 

v 7t 



which is the standard ground state of the oscillator obeying the vacuum condition 

aip (x, t ) = 0. 
The solution ipo(x, t) has the form of a Gaussian wave-packet 



ipo(x,t) 



1 1 



e 2=1*5 



(33) 



(34) 



(35) 



The Fock states ip n (x, t), which are solutions of the Schrodinger Eq.d3TI). are constructed by starting from the vacuum state 
by standard algebraic formulas using the integrals of motion ( f24b and (f25l 



Solutions d36*l > have the explicit form 



1p n (x,t) = 1po(x,t) 



m 

V2™^I \e(t) 



n/2 



H„ 



Mt)\< 



(36) 
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Here H n are Hermite's polynomials. There exist Gaussian packets which are squeezed coherent states which are obtained by 
means of the Weyl displacement operator acting on the vacuum state, i.e. 



-0q(x, t) = V(a)ipo(x, t) . 



where the Weyl system reads 



T>{a) = exp - a*A(t) ) 



Here a are complex numbers a = a\ + ia 2 and the squeezed coherent states have the properties 

|2 |fl|2 



J ^ a ( x , t )^(x,t)dx = exp(-^--^+a*p\ , 



and 



4>u{x,t)ip a (x' ,t)daia 2 = S(x - x') 



(38) 



(39) 



(40) 
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Let us now consider the tomographic probability description of the quantum parametric oscillator. The symplectic tomogram of 
the oscillator quantum vacuum state ipo(x, t) is expressed in terms of the wave function 



WoO,/x,^ t) 



1 



2tt \v\ 



Using the explicit expression of the wave function, we get a Gaussian tomographic probability distribution of the form 



W a (x,fi, v, t) 



1 



: exp 



X 2 
'2< 



where the dispersion of the random position X depends on the parameters fj, and v and the function e(t) as follows 



where 



\m\' 



* 2 pq = \{WW)\ 2 



(42) 



(43) 



(44) 



(45) 
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The state corresponding to the tomogram (1431 is the squeezed vacuum state. Depending on e(t) and e(t) the fluctuations of 
position or momentum can be either smaller than 1/2 or larger than 1/2. The state has the position-momentum correlation 
a pq 7^ and this correlation satisfies the minimization of the Schrodinger-Robertson uncertainty relation ll32l.l33ll 
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The coherent state ip a (x, t) has also a Gaussian tomogram 



W a (X,n,v, t) 



: exp 



2ct 2 



where the dispersion a 2 is given by Eqs. (144-b and (PES! and the mean value reads 



where 



Here 



X = n(q) a + v{p) a , 
(q) a = V2 Re a(t) , (p) Q = V2 Im a(t) . 



a(t) 
«*(*) 



V2 



V2 



e(t)p-e(t)q 
(e(t)p-e(t)i 



and 



V2Re 



a . 



p = \/2Im a , 



where a = a.\ + ia-i, i.e. a is a constant complex number labelling coherent states (x\a, t) = ip a (x, t). 
The probability distributions Wq{X 7 fi, u, t) and W a (X, fi, v, t) satisfy the kinetic equation 

VV(X, n, u, t) - ^W{X, /i, 1/, t) + w 2 (t)z/— n, v, t) = . 

The kinetic equation corresponds to the classical Liouville equation in classical mechanics, that is 

df(q,P,t) , df(q,p,t) df(q,p,t) dV(q,t) 



+ 



dt ' dq P dp 
where the potential energy V(q, t) is the energy for the parametric oscillator 



= 0. 



V{q,t) 



1 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



(54) 



However, the same Eq.(l52b corresponds to the quantum von Neumann equation for the density operator of the parametric 
oscillator 



being H = 1, where 



In the tomographic probability representation Eqs 
expressed in terms of the propagator 



i§jPm = [*(*). />(!)]. 



(55) 



(56) 



and (1551 1 coincide. The evolution of the tomogram W(X, fj,, v, t) can be 
W(X, n, u,t)= I U(X, fi, v, X\ n', 1/, t)W(X', fi', v 1 , 0)dX'dn'dv' . (57) 
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The propagator H(X, fi, v, X', //, v' , t) in this integral relation is identical for both classical and quantum evolutions of the one- 
mode-parametric oscillator with a time-dependent frequency. The difference appears in the initial conditions. However, classical 
tomograms have not to respect the uncertainty relation 



X 2 W(X,fi,v, t)dX\^ =1 , v=0 



XW(X,n,y, t) 



S^J X 2 W(X,fx,u,t)dX\^ v=1 



XW(X,^,u, t) 



//=(),;/=! 



1 



(58) 



while the tomogram of a quantum oscillator must satisfy these inequalities. This results means that, since the evolution is 
governed by the same equations, both situations, classical and quantum, are represented. Since the field is the collection of 
modes, the arguments presented above are applicable to all the fields. This means that for the field tomogram, which is the 
product of tomograms of all modes, the evolution result classical. Quantumness is hidden in the initial state but summing up all 
states gives, as result, the classical field. The perturbations described by the Hamiltonian, thus in the tomographic picture, can 
be discussed both using the quantum or the classical propagator of the field since both coincide. 

It is straightforward that the Hamiltonian describing the cosmological perturbations is the Hamiltonian of field harmonic 
vibrations. It is the natural Hamiltonian of a system whose properties fluctuate around the background state (the FLRW back- 
ground) and where the fluctuations are small (small vibrations). The deviations from the "equilibrium" background states are 
obeying the Hooke law. Only for large deviations from the equilibrium state, the Hooke law is violated and strong anharmonicity 
appears. The "spring" providing the harmonic vibrations of a system with very large amplitude of vibrations can even break 
giving possibility to the system to come out far away. This simple mechanism in the classical description of small vibrations and 
in quantum domain can give different results. In fact the "ground" state of the field is the squeezed state of a harmonic linear 
oscillators with time-dependent frequencies. The fluctuations of these vibrations in the classical domain, at zero temperature, 
are also equal to zero. The small vibrations cannot provide, in this case, large deviations from the equilibrium state even if 
one takes into account the possible appearance of Hooke's law violation for large amplitudes of the oscillations. For a quan- 
tum domain, even these small vibrations have quantum fluctuations. Thus the deviations of the oscillation position |e(i)| 2 /2 
and of the momentum |e(i)| 2 /2 can become large in comparison with the standard vacuum ones and be equal to 1/2, due to 
the influence of the possible large values of the contributions of |e(f)| (or |e(£)|). This mechanism of quantum fluctuations in 
presence of time-dependent frequencies can, for some parts of the system, create large amplitudes which violate the Hooke law. 
In such a case, the system can give rise to domains of larger densities and non-uniformities which results as inhomogeneities 
and anisotropies with respect to the background. This could be a coherent scheme to match quantum microscopic primordial 
perturbations with the today observed large scale structure. 



VI. THE MULTI-MODE PARAMETRIC SMALL PERTURBATIONS 



Let us consider now the most general multi-mode Hamiltonian. This situation results more realistic in order to take into 
account the problem of cosmological perturbations. As in yfl, we consider periodic conditions. Due to these, the Hamiltonian 
can be taken in the following form 

W(t) = \~QB(t)Q + C(t)Q . (59) 
Here we consider N modes (N can be equal to oo) and the vector operator 

Q = (Pi,P2, ■ ■ ■ ,PN,qi,q2, ■ ■ -,Qn) ■ (60) 

The system properties are coded by the interaction 2N x 2N matrix B(i) which depends on time t,that is 

The N x N block matrices -Bfe(t), k = 1, 2, 3, 4 correspond to quadratures interactions, namely Bi(t) describes i\— quadratures 
interactions, B±(t) corresponds to ^—quadratures interactions and the matrices B 2 (t) and B 3 (t) to interaction terms due to q 
and P couplings. The 2 TV-vector C(t) provides the interaction terms corresponding to homogeneous "electric-like" fields acting 
on charged particles. 
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The classical counterpart of this Hamiltonian, in the same quadratic form, is given by a vector ( f60b composed by classical 
momenta and positions with standard Poisson brackets. The solution of the Schrodinger equation with first integrals of motion, 
in position and momenta, determines an inhomogeneous symplectic group element parameterized by 2N x 2N matrix A(t) and 
a 2N— vector A(i). One has the 2N— vector I(t) whose components are integrals of motion. The vector reads as 



I{t) = A(t)Q + A(t). 



The symplectic matrix A(i) satisfies the equation 



(62) 



A(i) = A(i)£B(f), 



S = 



-liv 
Ijv 



(63) 



and the initial condition 



A(0) 



*- 2 N ■ 



(64) 



The vector A(t) satisfies the evolution equation 



A = A(t)EC(t) 



(65) 



with initial value, A(0) = . 

The evolution of the system of field modes with Hamiltonian (|59l is described by the Green function of the Schrodinger 
non-stationary equation determined in terms of the matrix A(t), the vector A, and the Gaussian Green function [24]. It is 



G{x,x',t) = [rfet[-27r. i A 3 ]- 1 / 2 



exp 



A 3 1 Xix — 2x\r, 1 x' + x'XiX 



+2x'Xz 1 5 2 + 2x(Si - Xi\£%) + % -2 \T%(r)dT 



(66) 



Here, the symplectic matrix Kit) as well as the vector A are in a block form being 



A(t) 



Xx{t) X 2 (t) 
A 3 (t) A 4 (i) 



A(t) 



5i(t) 



This means that, given the initial quantum state at t = 0, ip(x, 0) where x = (x\, x 2 , . . . xn), the state at time t reads 

ip(x,t)= / G(x , f, t)ip(x, t)dx' . 



If the initial state is given by a density matrix p(x, x' , 0), the density matrix at time t reads 

p(x,x',t) = f G(x,y,t)p(y,tf,0)G*{x',y',t)dydy' . 



(67) 



(68) 



(69) 

In ID], some initial states are studied. The first is that the initial state is the vacuum state with a Gaussian wave function 
corresponding to a set of oscillators with different frequencies. Another possibility is that the initial state is chosen as a thermal 
density matrix of the set of oscillators at temperature T. The density matrix is also Gaussian and the state at t = can be 
chosen as a pure squeezed state with a large, but finite, number of particles with wave functions expressed in terms of Hermite 
polynomials. The true vacuum state and the thermal state can be considered in both classical and quantum universe evolution 
pictures. The squeezed Fock states of the evolving Universe correspond to the quantum domain only because the states cannot be 
realized as initial state in the classical domain. Below we focus on the vacuum initial state to show, in tomographic probability 
representation, its evolution both in classical and in quantum pictures. 

The quantum tomogram W(X, (2, v, t), where X = (Xi, X 2 , . . . , X^), fl = (/ii, p 2 , . . . , /Ltjv), v = {v\, U2, . . ■ , vn) which 
is a joint probability density of random variables X/., k = 1, 2, . . . , N is expressed in terms of the wave function 



w{x,{i,p,t) 



i 



nf=i27rk| 



tp(y,t)exp < i 



' N 

E 

.*:=i 



2v k Vk 



XkUk 



Vk 



dy 



(70) 
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The Wigner function of the state is expressed as the wave function 



W(q,p, t) = / i> [q + $,t U* [q- ^t) e~^du. 



The tomogram ( fTOl i is connected with the above Wigner function by the Radon transform 



N 

W{X,n,v,t)= / W(q,p,t)H 
•* fe=i 



5 (X k - [ikqu- VkPk) 



dqkdpk 
2tt 



(71) 



(72) 



The density matrix in position representation satisfies the von Neumann equation 

.dp(x, x', t) 



at 



Here the operators Qg and Q~, are given by d60l ) where Pks 



\ (QsB{t)Q s ^j p(x, x', t) + C(t)Q s p(x, x', t) 
-\ (QjB(t)Q x ^ P(x, x', t) - C{t)Q-, p(x, x>, t) . 



d 

c):r,. 



Qkx -> x k and P k ~, 



The equation for the Wigner function (the Moyal equation, see [31:]) is given by ( f73l with the replacements 



(73) 

%£' ~* x 'k respectively. 



d 1 d 

dx k 2 dq k ox' k 2 dq k 

i d 



d 1 d_ 

% a 

2 2 <9p fc 



(74) 



The evolution equation for the tomogram ( TTOb or (l72l is obtained from the equation for the Wigner function by the replace- 
ments 



W(q,p,t) -► W(X,/x,z?,i) 



9 



(9 



d d 

V- 



dq k "dX k dP k dX k 



Ik 



d ( d 



dpk \dXk 



Pk 



a f o 



dv k \dX k 



(75) 



It is important to point out that the equation for the Wigner function and for the tomogram will coincide for both domains 
quantum and classical with the equations for probability densities f(q, p, t) satisfying the Liouville equation and for the quantum 
and the classical tomogram 



f n 

W(X, p, v, t) = / f(q,p, t) Y[ [S (X k - Pkqk ~ v k pk) dq k dp k ] , 
J fc=i 



(76) 



respectively. 

The tomogram is the probability density. Due to this one can introduce the Shannon entropy lfT6ll associated with the proba- 
bility density. The formula for the entropy reads 



Hi(0) 



W(X,p,P)lnW(X,p,if)doS. 



(77) 



p k = Sk cos k , v k = s fc 1 sin 9 k . 

It means that for evolution of the Wigner function W(q,p, t), in the quantum domain, and the probability distribution f(q,p, t), in 
the classical domain, the propagators are identical as well as propagators for tomograms both in quantum and classical domains. 
It is a property of systems with Hamiltonians of the form (|59l l and it corresponds to the Ehrenfest theorem. 
Thus, given the evolution of the tomogram in two states, the evolution of entropy can be defined as 



H0,t) = Hi0,t) + H 2 (6,t) . 



(78) 
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where 

6=(6 1 ,6 2 ,...6 N ), (79) 

and 

H 2 (e) = Hi(e + ^\, 8 + ^ = {e k + ^}. (80) 

One can introduce the sum entropy 

H(8,t)=H(8,t) + H^8 + ^,?j , (81) 

where 

= + (82) 

The entropy H (8, t) satisfies the inequality, see e.g. lfl6ll . 

H{9,t)> N\mre. (83) 
For the initial ground state of the universe, one has the saturation of the above inequality, i.e. 

H(8)=N\nne. (84) 

In the quantum domain, the entropy H(9) cannot be less than the value iVln7re. In the classical domain, the entropy H(6) 
characterizes the order or disorder in the field state. In the process of evolution of the Universe, this entropy is changing. 
The tomograms for classical states are probability distributions evolving for the quadratic Hamiltonian (|59l identically. If one 
measures the tomogram at the late time i, corresponding to the classical epoch, one can learn what was the initial state since at 
this period of time the tomogram corresponding to the classical state provides the possibility to calculate the probability density 
f(q, p, t), using the inverse Radon transform 

i r N 

f(q,P,t) = j-^ J W{X,^v,ty^=^-^»-^ Y[dX k d» k dv k . (85) 



In principle, this density can be used to describe the distribution of matter in galaxies 13411 and in clusters of galaxies 113 511 . 
Besides, one can introduce another characterization of the Universe state in terms of the tomogram. 



VII. CONCLUSIONS 



To conclude, we point out the main results of the paper. We have considered classical and quantum perturbations under the 
same tomographic standard. The small vibration Hamiltonian with time-dependent parameters can give account of the basic 
model which can be adopted in the quantum and in the classical regime. Specifically, the classical and quantum descriptions 
can be associated with the tomographic probability distributions. The classical initial states (classical tomograms) and quantum 
initial states (quantum tomograms) have different properties because the classical tomographic entropies can violate the quantum 
bound of the tomographic entropy. The tomographic entropy can be estimated at present epoch by measuring the space and 
momentum distributions of the matter in the universe. Since the propagator for both classical and quantum tomograms of the 
universe states is the same, one can trace back, in principle, the entropy data to the initial state of the universe. The observed 
data, extrapolated to the initial state by means of the classical tomogram, can be less than the quantum bound. The existence of 
bound property can be used to discriminate between cosmological quantum and classical behaviors. 

The physical meaning of the entropy which we have considered is that it represents the sum of two entropies. One corresponds 
to the spatial probability density of matter. Then, we have to assume 8 = 0. The other one corresponds to the momentum 
distribution. For this case, we have to assume the value 8 = ^. By measuring these distributions at present epoch, one can trace 
back the cosmological evolution to entropy initial values corresponding to the small vibrations. In other words, the tomographic 
description of the universe provides the possibility to correlate classical and quantum cosmological perturbations in the same 
unitary scheme. The final goal of this picture is to find out a dynamical and self-consistent approach capable of connecting 
the primordial quantum perturbations to the today observed large scale structure rf35tl . In a forthcoming paper, we will discuss 
the above results considering the data coming from the observations and their relations with initial values of cosmological 
parameters. 
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